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ABSTRACT 

We present observations of early-type galaxies NGC524 and NGC2549 with laser guide star 
adaptive optics (LGS AO) obtained at GEMINI North telescope using the NIFS integral field 
unit (IFU) in the /Tband. The purpose of these observations is to determine high spatial resolu- 
tion stellar kinematics within the nuclei of these galaxies and, in combination with previously 
obtained large scale observations with the SAURON IFU, to determine the masses (M.) of 
the supermassive black holes (SMBH). The targeted galaxies were chosen to have central light 
profiles showing a core (NGC524) and a cusp (NGC2549), to probe the feasibility of using 
the galaxy centre as the natural guide source required for LGS AO. We employ an innovative 
technique where the focus compensation due to the changing distance to the sodium layer is 
made 'open loop', allowing the extended galaxy nucleus to be used only for tip-tilt correc- 
tion. The data have spatial resolution of 0'.'23 and 0'.'17 FWHM, where at least ~ 40% of flux 
comes within CC'2, showing that high quality LGS AO observations of these objects are pos- 
sible. The achieved signal-to-noise ratio (S/N~ 50) is sufficiently high to reliably determine 
the shape of the line-of-sight velocity distribution. We construct axisymmetric three-integral 
dynamical models which are constrained with both the NIFS and SAURON data. The best 
fitting models yield M.=(8.3± 2 { 7 3 ) x 10 8 M Q and (M/L)i = 5.8 ± 0.4 for NGC524 and 

M.=(1.4±5 j) x 10 7 M and (M/L) R = 4.7 ± 0.2 for NGC2549 (all errors are at the 3a 
level). We demonstrate that the wide-field SAURON data play a crucial role in the M/L deter- 
mination increasing the accuracy of M/L by a factor of at least 5, and constraining the upper 
limits on black hole masses. The NIFS data are crucial in constraining the lower limits of M. 
and in combination with the large scale data reducing the uncertainty by a factor of 2 or more. 
We find that the orbital structure of NGC524 shows significant tangential anisotropy, while at 
larger radii both galaxies are consistent with having almost perfectly oblate velocity ellipsoids. 
Tangential anisotropy in NGC524 coincides with the size of SMBH sphere of influence and 
the core region in the light profile. This agrees with predictions from numerical simulations 
where core profiles are the result of SMBH binaries evacuating the centre nuclear regions fol- 
lowing a galaxy merger. However, being a disk dominated fast rotating galaxy, NGC524 has 
probably undergone through a more complex evolution. We test the accuracy to which M. can 
be measured using seeings obtained from typical LGS AO observations, and conclude that for 
a typical conditions and M. the expected uncertainty is of the order of 50%. 

Key words: galaxies: individual (NGC524, NGC2549)- galaxies: elliptical and lenticular - 
galaxies: kinematics and dynamics - integral field spectroscopy 



1 INTRODUCTION 

The discovery that the masses of central dark objects correlate with 
the global characte ristics of their host galaxies, such as the total 
lumin osity or mass jKormendv & Richstonell 19951 : iMagorrian et alj 
1998), led to a realisation that these central objects are cru- 
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cial for our interpretation of galaxy formation and evolution 
(e.g. ISilk & Reeslll998l; |Pi Matteo et atll2005l : ICroton et alj|2006l: 
iBower et al.l 120061 : iHopkins et alj|200q) . The mass of these dark 
objects, and their link with active galactic nuclei suggest they are 
supermassive black holes (SMBHs). The most secure evidence 
for the existence of a nuclear SMBH comes from the centre of 
the Milky Way, where the large mass is determined by the three- 
dimensional orbits of stars rotating around it (e.g. Schod eTet alj 
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l2002l : lGhez et al.l2003ll200^ : lGillessen et alj2009h . Its small phys- 
ical size is confirmed by the variation of spatially unresolved ac- 
tivity which is linked to an accretion disk around the obj ect (e.g. 
IDoeleman et alj200ll : Ishen et alj2005t iFalcke et al.ll2009h . 

Such detailed information from close to the black hole is, 
however, not available for more distant non-active galaxies which 
exhibit the scaling relations between the M. and the host galaxy 
properties. The masses of these SMBHs can be determined using 
dynamical tracers such as a disk of gas clouds in Keplerian ro- 
tation around the SMBH or the line-of-sight velocity distribution 
(LOSVD) of stars coupled with dynamical models, which account 
for the complexity of possible stellar distribution functions. Both 
methods require observations of high spatial resolution since the 
models have to be constrained by resolved kinematics which probe 
the sphere of gravitational influence of the SMBH (R sp a), denned 
as the distance from the black hole at which the potential of the 
galaxy and SMBH are approximately equal. This spatial scale is 
usually defined as R Bp h = GM./a 2 , where a is the velocity dis- 
persion of the galaxy (usually measured over a large aperture), and 
is typically significantly less than an arcsecond in apparent size, 
even for nearby galaxies. 

While K sp h is not a hard limit, and the influence of the 
SMBH will be felt to a lessening degree at larger radii, it does 
provide a useful estimate of the scale on which the observable 
effects of an SMBH on the galaxy's dynamics will be most sig- 
nificant. Until recently, only long-slit observations from Hubble 
Space Telescope (HST) were able to achieve this resolution. These 
observations greatly improved the accuracy of the M. determi- 
nations and led to a discovery of eve n tighter relations be tween 
M. and a e dFerrarese & Merrittll200(il ; iGebhardt et alj|2000h fol- 
lowed by similar correlations between M. and host galaxy prop 
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Even the HST observations, however, have their limits. The 
HST is a 2.4 m telescope, has only a long-slit spectrograph and it 
is optimised to work at optical wavelengths. These specifications 
prohibit the use of the telescope for measuring spatially resolved 
kinematics in small, faint galaxies, in low surface brightness giant 
galaxies with cor e-like nuclear p r ofiles and in galaxies w ith dust- 
obscured nuclei (Ferrarese 20031; iFerrarese & Fordll2005h . Still, a 
decade of observations with the HST resulted in about 45 secure M. 
determinations, mostly of early-type galaxies with a e = 150 — 350 
km/s and spanning th ree orders of magnitude in M. dGrahamll2008l : 
iGultekin et alj2009h . 

The advent of adaptive optics (AO) facilities at the 8- 10m 
class of ground based telescopes opens a new door into the study 
of SMBHs. These observations, whether natural guide star (NGS) 
or laser guide star (LGS) assisted, are approaching the resolution 
of the HST at near-infrared wavelengths which are also less af- 
fected by dust and, hence, can be used to probe dusty nuclei of 
jiral galaxies extending the demography of the M. — a relation 
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Cappellari et alj|2009l) . In addition, several of the instruments be- 



hind these AO systems are integral-field units (IFUs), which, given 
that they map the LOSVD over a two-dim ensional area, offer a 
tighter constraint on the orbital distribution dVerolmeet all [2002) 
needed for a robust recovery of distribut i on function from the 
observables dCappellari & McDermidll2005l ; iKrajnovic et alj|2~005l ; 
Ivan de Ven et alj|2008l) . 

The link between the SMBHs and the formation of the host, 



both in terms of mechanisms and the evolution of this link, is still 
somewhat hidden behind the uncertainties of the M. — a relations. 
The shape, extent and scatter of this relation depend not only on the 
accuracy of the M. determinations, but also on the usage of a con- 
sistent measure of a, and the selection biases of the galaxies used in 
this relation. In this respect, future efforts need to be focused on ex- 
tending the range of a probed in this relation, ideally for unbiased 
samples using consistent techniques for measuring M.. This work 
aims at taking the first step in this direction, refining our ground- 
based observational techniques to allow larger numbers of objects 
to be studied in a homogeneous way. 

Spectroscopic observations of nuclei can also be used to probe 
the formation mechanism of SMBHs. If central SMBHs are ubiq- 
uitous in galaxies (or at least present in those with bulges), then 
when smaller galaxies merge, the two central black holes sink to 
the bottom of the potential well forming a binary which should 
also merge over some time. This process can only be completed 
with the assistance of nearby stars; the energy and the angular mo- 
mentum of the binary is transferred to stars approaching the binary 
on radial orbits, which are then ejected at much higher vel ocities 
dOuinlan & Hernquistl 1 19971; IMilosavlievic & Merrittl 1200 ll) . This 
process leads to dynamical heating of the centr al regions and a 
decre ase in the nuclear surface brightness profile dKormendv et al.l 
2008). The 'scouring' of the core is used to explain the existence of 
core-like light profiles (Faber et al. 1997), a nd nuc l ear mass deficits 
in massive ellipticals (Milosavli evic et all |2002| ; iKormendv et al.l 
120081 ; IKormendv & Benderll2009l ). In addition to producing a core 
in the light profile, simulations of this process also predict that 
the merging of an SMBH binary would leave a clear signature 
in the orbital distribution: there sho uld be a bias towards tan- 
gential orbits in the core region (e.g. iMakino & EbisuzakHll996l ; 
lOuinlan & Hernquistlll997l; IMilosavlievic & Merrittl200ll) . 

The main prohibitive issue with AO observations is that, even 
with LGS capabilities, a natural guide star is required for an optimal 
correction of atmospheric aberrations. There are only a handful of 
galaxies with a near-resolvable R sp h that have a suitable guide star, 
even at the fainter magnitudes permitted by LGS AO. Moreover, the 
correction provided by current single reference source, single con- 
jugate AO systems quickly degrades when the off-axis distance of 
the guide star exceeds the isoplanatic patch size, which is typically 
much smaller than the effective radius of the galaxy. This means 
that the use of stars for AO observations of nearby galaxies is not 
practical for large samples. 

This limitation has prompted innovative use of AO techniques, 
such as neglecting altogether the low-o rder corrections provided by 
the natural tip-tilt dDavies et al. 2008). In this study, we employ a 
novel method developed at Gemini Observatory which allows the 
galaxy centre itself to be used as a reference for tip-tilt correction. 
With this study, we explore the performance of this technique using 
two objects with differing central light profiles: NGC524 which has 
a flattened 'core' profile, and NGC2549 which has a steep 'cusp' 
profile. These two representative cases will demonstrate the feasi- 
bility of applying this technique to larger samples of objects, which 
in turn will allow us to tackle the issue of refining and better under- 
standing the connection between SMBHs and their host galaxies. 

In Section|2]we present the observation, in particular the LGS 
setup and the data used. In Section [3] we describe the data reduc- 
tion and extraction of the kinematics. In Section [4] we discuss the 
determination of the PSF and AO correction, while in Section[5]we 
present the dynamical models and the main results of this paper. In 
Section[6]we discuss and in Section|7]we summarise our results. 
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2 OBSERVATIONS AND LGS AO SETUP 

The results of this work are based on the combination of large- 
scale integral-fiel d kinematics from the SAURON instrument 
de-acon et al.ll200lh and high spatial resolution integral-field kine- 
matics from NI FS on Gemini North. The SAURON observations 
are presented in lEmsellem et al.ld2004l) . In this work we present the 
new data obtained with NIFS on Gemini North using laser guide 
star adaptive optics (LGS AO). 

2.1 LGS AO observations 

The Altair LGS system consists of a laser that projects a ~ 10 Watt 
coherent beam at 589 nm into the sky, exciting the sodium atoms at 
the height of around 90 km in the atmosphere, which then re-emit 
radiation and glow as a source on which Altair can perform the 
high-order corrections. A more complete d escription of the Altair 
LGS system is given in Bo ccas et ai]( l2006b . 

The final improvement of the PSF with the LGS system still 
depends on the existence of a nearby (< 25" for Altair) natural 
guide star (NGS). This star is needed to correct for the angle-of- 
arrival variations due to the atmospheric turbulence, usually called 
the tip-tilt correction, and to measure the absolute focus of the tele- 
scope. The first correction arises because the deflection of the beam 
of the LGS when it travels upwards and downwards through the 
atmosphere cancel each other, while the beam of the object above 
the atmosphere is deflected only once, leaving the low-order, tip- 
tilt, disturbance not detected by the LGS AO system. The second 
correction is a consequence of the finite and slowly varying dis- 
tance between the telescope and the sodium layer, which prevents 
the use of the LGS beacon for determining the focus of the tele- 
scope at infinity. Both degradations, tip-tilt and change in focus, 
can be corrected with an NGS (often called 'tip-tilt' star), which 
can, in general, be further away and a few magnitudes dimmer (up 
to 3 magnitudes for the Altair LGS system) than for a pure NGS 
AO observation. 

For a given guide-source flux rate, LGS AO typically yields 
smaller Strehl ratios than NGS AO observations due to the effects 
of the finite size of the LGS beacon in the atmosphere, the cone 
effect and increased sensitivity to the actual seeing, which distorts 
the upward travelling laser beam. For bright, on-axis sources, the 
predicted Strehl ratio can b e as much as 20% les s with an LGS 
system compared to NGS £e Louarn et al. I Gil). Nevertheless, 
they dramatically extend the useful sky coverage, allowing AO- 
impro ved observations of targets for which a bright NGS is unavail- 
able dStuik et al.l200j) . The Altair systems requires, for a low Strehl 
observations, a tip-tilt star of 17.5 (R band) magnitude to be within 
25"from the target, the nucleus of the galaxy in our case. Very few 
galaxies with R sp h > 0'.'05 (typically nearby early-type galaxies) 
have neighbouring stars on the sky that meet this criterion. Increas- 
ing the accessible off-axis distance only degrades the PSF, and the 
use of fainter stars becomes problematic due to the non-negligible 
background of the extended galaxy at these magnitudes. This poses 
a serious limitation on what can be achieved on this topic with ex- 
isting single-guide source, single conjugate AO systems on large 
telescopes. 

2.2 The open-loop focus model 

To avoid the issue of suitable guide stars, it is, in principle, also 
possible to use the nucleus of the galaxy as a natural guide source 
for the system, provided it is sufficiently bright and compact. For 



Altair, the 'rule of thumb' is that the nucleus should show a drop 
of > 1 magnitude in brightness within the central ~ 1", which 
ensures a certain degree of contrast within the field of view of 
the wave-front sensor. This condition is met for many AGNs or 
qu asars, which have bright, unres olved nuclei (e.g. the case of Cen 
A. lHaring-Neumaver et alj|2006l) , and many nearby spirals have a 
suitably bright and distinct central star cluster to allow AO correc- 
tions. But for quiescent early-type galaxies, the situation is less ob- 
vious. The central light profiles of early-type galaxies show a gen- 
eral change from steeply rising 'cusp' profiles at lower masses, to- 
wards a flatter central profile that can de fine a central 'core' region 
dFerrareseetai]|l994 iFaber et alj| 19971) . For the most massive and 
nearby quiescent galaxies, the core region (inside which the surface 
brightness changes only little) can be larger than the wave-front 
sensor field of view, which sees essentially a constant background 
of light from which no tip-tilt information can be obtained. How- 
ever, this is only relevant for galaxies at the most massive end of 
the M. — a relation. The majority of galaxies have either core radii 
that are similar to or less than one arcsecond, or exhibit steep, cuspy 
profiles which should show contrast within the AO wave-front sen- 
sor. In order to test what correction of the PSF could be expected 
for typical early-type galaxies, we chose NGC524 and NGC2549 
which both satisfy the basic Altair requirement for n uclear guid- 
ing. Specifically, NGC524 has a core-like light profile dFaber et alj 
1 19971). while NG C2549 a cusp-like profile at the HST resolution 
dRest et al .11200 ll) . 

After initial observations were attempted, it became clear that 
the chosen nuclei, while suitable for tip-tilt correction, are too faint 
to be used for constraining the focus, resulting in an unstable fo- 
cus control loop and subsequently compromised image quality. To 
avoid this, Gemini staff implemented a procedure by which the fo- 
cus correction during the science integration is controlled by a geo- 
metric function that takes into account the change in the distance to 
the sodium layer as the telescope position changes. This, so-called 
'open-loop' focus model, was a pre-existing feature of the LGS- 
AO system, to reduce the convergence time of the system following 
large slews during night-time operations. The only time-dependent 
parameter of this model is the altitude (i.e. absolute height above sea 
level) of the sodium layer, which is determined immediately before 
observing the galaxy by 'tuning' the LGS AO system using a nearby 
bright star. When all the control loops of the system (tip-tilt, focus, 
and LGS) have converged with this reference source, the loops are 
opened and the science target is acquired. The tip-tilt and LGS con- 
trol loops are then closed, using the galaxy nucleus and laser beacon 
respectively as reference sources. The focus loop is left open, being 
passively controlled by the open-loop model. After approximately 
one hour of observations of the galaxy nucleus, the bright reference 
star is re-observed; first with the focus loop left open, and then with 
the focus loop closed, in order to measure the relative degradation 
of the PSF due to any departure from the open-loop model during 
the science integration. Consequently, the LGS AO system is now 
re-optimised at this stage for further observations. 

We note that, whilst this technique was under development at 
the time our observations were obtained, this mode of operation is 
now fully supported, as of the 2009A semester. 

Figure[T]shows the stability of the PSF correction in the open- 
loop focus model. The panels show horizontal cuts through the re- 
constructed image of the tuning star observed with LGS AO with an 
actively optimised focus (left panel) and taken one hour later, where 
the focus of the telescope has been completely controlled during 
that time by the open-loop model (right panel). The observed degra- 
dation of the PSF is minimal. In both cases, the measured PSFs 
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Figure 1. The stability of the PSF correction in the open-loop focus model 
as measured on the horizontal cuts through the reconstructed images of a 
star used to set up the LGS AO system. On both panels blue and orange 
lines are the narrow and broad Gaussian components describing the PSF. 
The red line is the sum of these two components. Left: The LGS AO is 
fully tuned in, with both tip-tilt and focus loops closed and using the star 
for guiding. The FWHMs are 0.' 1 1 and 0.'48, for narrow and broad com- 
ponents, respectively. Right: One hour later, the LGS AO system is in the 
same set up as during the observations of the science target: tip-tilt loop is 
closed on the star, while the focus loop is open and the changes in the focus 
still follow a model. The degradation is virtually negligible: FWHMs are 
0.' 12 and 0.'60 for narrow and broad components, respectively. 

can be well approximated by a double Gaussian, with full-width 
half maximum (FWHM) values changing from O'.'ll and 0'.'48 to 
f/.'l2 and 0'.'60 for the narrow and broad components of the PSF, 
respectively. Approximating the Strehl ratio as the ratio of the to- 
tal flux in the narrow Gaussian component and the total flux in the 
whole PSF, we obtain the values of ~ 31% for both closed and 
open loop PSFs. These negligible differences show that it is possi- 
ble to achieve a very good PSF correction of 0" 1-0.' 2, (FWHM) ap- 
proaching the diffraction limit in the K band, even without the focus 
loop closed. On this occasion at least, the Hawaiian skies were sta- 
ble over approximately one hour, demonstrating that the open-loop 
focus model used by the Gemini system seems to track accurately 
the changes in distance to the sodium layer. 

The PSFs measured from the tuning star observations, how- 
ever, cannot be considered representative of the PSF of science ob- 
servations obtained guiding on the nucleus. The above demonstra- 
tion reassures us that the change in focus during the galaxy observa- 
tion was well approximated by the open loop model, and suggests 
that the delivered PSF of the galaxy centres may also comprise a 
narrow component comparable to the diffraction limit and a sub- 
stantial broad component. Estimation of the actual delivered PSF 
for our galaxies, with their differing central light profiles, is dis- 
cussed in Section[4] 



2.3 NIFS observations 

NGC524 and NGC2549 were observed during the 2007B and 
2008A semesters in the K band (central wavelength is 2.2/^m) with 
H+K filter and spectral resolution of R~5000. The field of view 
of NIFS is ~ 3 x 3". The observations were divided into separate 
blocks of approximately one hour with 600 seconds science expo- 
sures set in the sequence: OSOOSO, where O is an observation of 
the object (galaxy) and S is an observations of the empty sky field. 
NGC524 and NGC2549 were observed in total for 21 x 600 and 
15 x 600 seconds (3.5 and 2.5 hours) on source, respectively. Some 
frames had to be, however, discarded in the final stage because of 
inferior PSF corrections (see Section [4] for more details) and the 



final data cubes for NGG524 and NGC2549 were constructed out 
of 10 x 600 and 11 x 600 on source observations, respectively. In 
addition to galaxies, we observed telluric standards of A2 V and 
GO V types, switching between the stars and the sky fields in the 
same sequence as for the science targets (see Section [3] for telluric 
correction). 

2.4 Other data 

In the construction of dynamical models next to kinematic data 
we also use images of our galaxies. Since we need to constrain 
the models at the very high spatial resolution in the centre and 
at large radii in order to collect all available information about 
the distribution of light, we make use of Hubble Space Telescope 
(HST) images obtained with Wide Field and Planetary Camera 2 
(WFPC2) in F702W (NGC2549) and F814W (NGC524) bands, and 
ground based images obtained at the MDM telescope in F814W 
band (Falcon-Barroso et al., in prep). 



3 DATA REDUCTION AND STELLAR KINEMATICS 
EXTRACTION 

NIFS data were reduced using the set of reduction routines pro- 
vided by the GEMINI and incorporated within the IRAF package, 
following nifsexamples scripts for the calibration, telluric and sci- 
ence observationo Our data reduction departed from the GEM- 
INI guidelines in two instances: at the preparation of the telluric 
stars and co-addition (merging) of the individual fully reduced (flat 
fielded, sky subtracted, wavelength calibrated, telluric corrected) 
science frames. In this section we first describe these two steps and 
afterwards we outline the method we used for the extraction of kine- 
matics. 

3.1 Telluric correction 

Near-infrared ground-based spectroscopic observations are charac- 
terised by the presence of strong atmospheric absorption features. 
These features exhibit temporal and spatial dependency and have to 
be removed from the spectra in order to uncover intrinsic spectral 
characteristics of observed objects. The correction method is based 
on the fact that by observing a featureless star, both near in time 
and sky coordinates to the object one studies, one can obtain a ref- 
erence spectrum of the unwanted absorption lines. While seemingly 
straightforward, this method is hindered by the lack of objects that 
are both featureless, reasonably bright and evenly distributed on the 
sky such that can be found sufficiently near (in air mass) to poten- 
tial targets. Usually, stars of type A are used since they generally do 
not have strong metal features, show little reddening and are well 
appr oximated by a black b ody spectrum of 10 000 K. 

iMaiolino et alj < fl996h outlined how F and G stars can also be 
used if a high resolution solar spectrum is employed to remove the 
intrinsic features. In this case a telluric correction spectrum is de- 
rived dividing the observed G or F star by a solar template spectrum 
shifted and broadened to the appropriate relative velocity and res- 
olution. The quality of the correction depen ds on how similar the 
telluric star is to the sun. IVacca etHEl d2003[) put forward a simi- 
lar method that is based on AO V stars, such as Vega. These are a 
few times more common the G2 V stars and, hence, more likely 

1 http://gemini.conicyt.cl/sciops/instruments/nifs/NIFSReduction.html 
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to be found next to science targets, and can be used to derive a 
telluric correction spectrum by dividing them by a high-resolution 
normalised model spectrum of Vega, again shifted and convolved to 
the resolution of the telluric stars. 

We followed these two methods and observed two types of tel- 
luric stars: a GO V Hip3033 during the observations of NGC524 and 
A2 V Hip44717 and Hip31665 stars for NGC2549. Hip31665 was 
observed only one night (March 29 2008) covering 4 exposures of 
NGC2549, the rest begin covered by Hip44717. The telluric stars 
were reduced as science frames (sky subtracted, fiat-fielded, spa- 
tially and spectrally rectified) and one-dimensional spectra were ex- 
tracted . As in the above studies, we use a high-resolution solar spec- 
trurr0 ( Livingston & Wall ace! 1 99 lh and a similarly high-resolution 
model of Vegqj as templates to establish the relative velocity shift 
and broadening of the telluric st ars. We do this using the penal ised 
Pixel Fitting (pPXF) method of ICappellari & Emselleml d2004 see 
Section [3~3"1 for more details about the method), which conveniently 
allows us to select regions of the observed spectrum which are rel- 
atively free of telluric absorption, and derive the velocity shift and 
broadening (described by a truncated Gauss-Hermite polynomial) 
required to match the star's intrinsic absorption features. Since the 
high-resolution template spectra are flux calibrated (the normalised 
solar reference spectrum was renormalized using a black-body tem- 
perature of 5800 K), no polynomial term is included in the pPXF, 
so the residuals of the fits provide both the telluric correction and 
flux calibration spectrum. We then returned to the GEMINI NIFS 
pipeline and used the derived telluric correction spectra with task 
NFTELLURIC to correct each individual data cube for telluric ab- 
sorption features. 



3.2 Merging of science frames and binning of spectra 

Before individual science frames can be combined into a final data 
cube it is necessary to put them on a common wavelength range 
and sampling, and determine the relative positions (or to re-centre 
the frames). Frames with bad PSF should also be excluded from the 
combined data cube, and we discuss this in Section|4] We re-centred 
the science frames for each galaxy in the following way: we choose 
a frame with a reasonable PSF (judged by eye to be one of the 
best and showing circular and concentric isophotes in the very cen- 
tre), and compared all other frames to it. The comparison was done 
on the isophotes of the reconstructed images (created by summing 
in the spectral direction), such that the isophotes of the re-centred 
frames overlap, where the goodness of the overlap is judged by eye. 
Given the variation in the PSF correction, is it not always possi- 
ble to achieve a perfect match between the isophotes, especially the 
central ones, where the influence of the seeing is most significant. 
Centering on the outer isophotes is usually robust, and the typical 
uncertainty in re-centering was around f/.'Ol. After all frames were 
re-centred, we merged them into a final data cube. This was done 
by a custom-made IDL procedure that reads in all individual (re- 
cantered) frames, determines their spatial extent, defines a new gird 
of a given pixel size, and interpolates individual frames to this new 
common grid. The flux values of the final data cube were obtained 
as the median of the individual cubes. For both galaxies we interpo- 
lated to a grid covering 3" x 3" with 0'.'05 square pixels. These pix- 
els over-sample in the cross-slice direction, which is justified by the 
oversampling provided by our dither pattern, which uses step sizes 



equivalent to non-integer numbers of slices, and the large number 
of individual exposures included in the merge. We slightly reduce 
the sampling along the slices, although the PSF is still adequately 
sampled. 

Extraction of the higher-order moments of LOSVD requires 
signal-to-noise ratio, S/N, of generally high values, e.g. S/N > 4 
( Ivan der Marel & Franxl 1 19931: iBender et alj Il994l; IStatlerl [l995h . 
We use the Voroni binning technique o f lCappellari & Copinl d2003l) 
to ensure a minimum S/N whilst optimising the spatial resolution. 
We first estimate the noise of the unbinned spectra as the standard 
deviation of the difference between each spectrum and its median 
smoothed version (over 30 pixels). For both galaxies we require a 
target S/N of 80, which results in 370 and 377 bins for NGC524 
and NGC2549, respectively. In the central 1 arcsecond of the NIFS 
field-of-view, the data are binned to no larger than O'.'l x O'.'l to- 
tal size. We checked a posteriori the achieved S/N in each bin by 
comparing the median value of a spectrum with the standard devia- 
tion of the residuals (difference between the optimal stellar template 
and the observed spectra, see Section [33}. In most cases for both 
galaxies, the final S/N per bin was in the range from 40-60. 



3.3 Stellar kinematics 

We use the penaliz e d Pix el Fitting (pPXF) method of 
Cappellar T& Emselleml {2004) to derive the LOSVD in pixel 
space , parameterised by a Gauss- Hermite polynomial {Gerhard 
1 19931 : Ivan der Marel & Franxl fl993h . pPXF finds the best fit to a 
spectrum by convolving a template spectrum with the correspond- 
ing LOSVD given by the mean velocity V, velocity dispersion a 
and Gauss-Hermite moments hz and higher, which can be used to 
quantify the asymmetric and symmetric deviations of the LOSVD 
from a Gaussian. 

To avoid template mismatch effects, an 'optimal template' 
spectrum is derived as the linear combination of spectra from a tem- 
plate library which, for a given LOSVD, best represent the galaxy 
spectrum. Ideally, the template library should contain stars which 
are representative of the galaxy's stellar population. The K band 
light of early-type galaxies is dominated by cool and evolved gi- 
ant stars (red gia nt and asymptotic giant branch stars). We used the 
stellar library of Wing e et alj ( 120081 version 1.0, April 2007) pub- 
licly available from the GEMINI NIR Resources web sitfl This 
library consists of 29 stars covering the spectral types from F7 to 
M3 and luminosity classes I, II, III and V. We use a subset of these 
stars that was observed in the (combined) wavelength range from 
2.15 to 2.43 /im. Although the stars were observed with GNIRS, 
the library contains reduced one dimensional spectra with instru- 
mental dispersion Oinst = 26 km/s, which compares well with the 
spectral resolution of our data: Oinst = 30 km/s, as measured from 
the broadening of night-sky emission lines. The publicly available 
spectra are not corrected to rest wavelength, and must be brought to 
a common velocity before they can be combined as a velocity tem- 
plate. This was done using pPXF to measure relative shifts between 
the library stars, using the high-resolution solar spectrum used in 
Section l3~T1 as a template. The distribution of relative velocities was 
checked using a star from the library itself as a template, to explore 
the effects of template mismatch. After de-redshifting the spectra 
with the solar template, the dispersion in velocities measured using 
the library template was less then 2 km/s. 

In the K band, the stellar kinematics is mostly constrained by 



2 ftp://nsokp.nso.edu/pub/atlas/photatl/ 

3 Obtained from lKuruc3 jl991): http://kuracz.harvard.edu/stars.html 



http://www.gemini.edu/sciops/instruments/nifs/near-ir-resources 
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Table 1. Stellar templates used in the construction o f optimal templates for 
NGC524 and NGC2549 from the template library of lwinee et all 120081) . 



2.15 2.20 2.25 2.30 2.35 2.40 2.45 

X M 

Figure 2. Optimal templates for NGC524 (top) and NGC2549 (bottom). In 
both cases, black solid lines show the total galaxy spectrum (sums of spec- 
tra of all spatial bins), red lines are the optimal templates, while green dots 
are the fit residuals, shifted upwards by an arbitrary amount. Regions with- 
out residuals were excluded from the fit (only for NGC2549). Vertical lines 
and associated names on the top panel describe the main absorption lines 
detectable on spectra of both galaxies. Note that the (5-3) 12 CO feature is 
fairly well reproduced even if it lies outside the fitting range for NGC2549. 



the CO band head starting at A = 2.29 /xm and consisting of CO 



(2,0), (3,1), (4,2) and (5,3) transitions ^Wallace & Hinkldil997t) . In 
our spectra the information on the stellar continuum can only be 
found blueward of the 12 C0 (2,0) transition, and while it is impor- 
tant to keep as much as possible of this region during the extrac- 
tion, the blue end of the usable spectra is given by the blueward 
extent of stellar templates. In this largely featureless region possi- 
ble absorption-li nes are Na I, Fe IA, Fe IB, Ca I and Mg I (e.g. 
ISilva et alj2008h . 

Before calling pPXF we rebin all spectra in wavelength to a 
linear scale x=ln A, while preserving the number of spectral pixels, 
and convolve the template spectra with the (quadratic) difference in 
the spectral resolution between the NIFS and GNIRS spectra. We 
describe the LOSVD with the V, a, h$ ~ ha moments, and add a 
4th order additive polynomial to correc t for the shape of t he con- 
tinuum, given that the template stars o f IWinge et alfd2008|) are al- 
ready continuum subtracted. Prior to extracting the kinematics, we 
performed Monte Carlo simulations to determine the level of penal- 
isation used, and determined that A = 0.4 is op timal for our data 
(see definition in ICappellari & Emsell em 2004). The penalisation 
reduces the scatter in the derived kinematic parameters when the 
intrinsic dispersion becomes similar to the instrumental dispersion 
(and therefore usually critically sampled) by biasing the LOSVD to- 
ward a simple Gaussian. In principle, this is more important in the 
case of NGC2549 t han for NGC524 since their a e are 145 and 235 
km/s, respectively dEmsellem et alj|2007l) . In both cases, however, 
the LOSVDs are expected to be well sampled by the observational 
set up (a e > pixel scale in km/s), making the penalisation less im- 
portant. 

Rather than determining an optimal template for every bin, a 



NGC524 



NGC2549 



star 


type 


star 


type 


_HX>113538t 


K8V 


ffX>113538 


K8V 


_ffD32440t 


K6III 


HD2490 


M0III 


HD63425B 


Kill I 


HD4730 


K3III 






HD63425B 


K7III 






HD720 


K5III 



Notes - Information on the type was obtained from GEMINI NIR Resources 
web site. Stars with f are given different type in the Simbad Astronom- 
ical Database I http://simbad.u-strasbg.fr/simbad/ 1: HD113538 is K9V and 
HD324440 is K4TO 



'global' optimal template was derived for each galaxy, which was 
then held fixed for each bin in that galaxy (the polynomial is still 
free to vary - only the relative mixture of spectra from the library 
is fixed). Each optimal template was determined from a pPXF fit to 
the sum of all the binned spectra of the galaxy, using the full tem- 
plate library of 23 stars. The resulting optimal templates show cer- 
tain differences in the stellar populations between the two galaxies. 
While for NGC524 we considered the full wavelength region com- 
mon to the stars and the NIFS spectra spanning 2.18 — 2.42/im, this 
approach for NGC2549 yielded a considerable template mismatch, 
mostly visible in the spatially non point-symmetric distribution of 
hz values across th e field-of-view - oft en taken as an indicator of 
template mismatch (Bender et al. 1994). In addition, (2-0) and (5- 
3) 12 CO absorption features were not reproduced well, their depths 
being under- and over-estimated, respectively. In order to minimise 
this template mismatch, we masked Nal and Cal absorptions lines 
and truncated the fitting region before the fourth 12 CO (5-3) tran- 
sition when fitting the combined spectrum of NGC2549. Figure [2] 
presents the final optimal templates, the absorption lines used and 
the residuals to the fits. 

The resulting optimal template for NGC524 combines only 
three stars from the library, while in the case of NGC2549, a to- 
tal of five stars were selected by the pPXF fit to form the optimal 
template. The stars, their types and the weights are presented in Ta- 
ble[3] In both cases, a K giant star can reproduce the main spectral 
features (the CO bandhead) to the first order, but to fit the full spec- 
tral range a cool dwarf (HD1 13538 in the library) is also needed. 
The dwarf star is particularly important for reproducing spectral 
features blue-ward from the CO bandhead, especially the promi- 
nent Na I (2.20pm), Ca I (2.26pm) and Mg I (2.28pm) absorptions. 
This star, however, does not reproduce well the CO absorption fea- 
tures. The remaining small discrepancies are visible mostly in the 
fit to the (2-0) transition, and in the case of NGC524 for the Nal 
absorption feature, while weak 13 CO features are reproduced well 
(although they are almost blended with the dominant CO bandhead 
due to the larger velocity dispersion, and are partially outside the 
fitting range for NGC2549). 

After the optimal template is constructed, we apply it to the 
spectra of the individual bins, using the same wavelength regions 
during the fit as in the construction of the optimal template. While 
it is possible that there are variation in stellar populations between 
bins which can not be described by the optimal template, the re- 
maining difference can be well fitted by the additive polynomials 
alone. We visually inspected all spectra to verify that they are well 
reproduced by the given optimal templates and we recheck if there 
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Figure 3. Kinematics maps of NGC524 (left) and NGC2549 (right). From 
top to bottom maps present: the mean velocity V, velocity dispersion <r, 
and Gauss-Hermite moments: /13, /14, /15 and hg. The colour-bar at the 
bottom show the plotted colour range, while the lower and upper limits for 
each map are given to the right of NGC2549 maps as set of two numbers. 
In the case of a maps, the upper pair of numbers refers to NGC524 and 
lower to NGC2549. Note that even /15 and he, although noisy, retain their 
main characteristics of being anti-symmetric and symmetric across maps, 
respectively. North is up and East to the left. 
© 2007 RAS, MNRAS 00Q,[THT9l 
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Figure 4. Comparison between NIFS (red open circles) and SAURON 
(filled circles) obtained by averaging data within concentric circular rings 
of different radii for NGC524 (top) and NGC2549 (bottom). NIFS data for 
NGC524 were systematically lowered by 4% to match the SAURON data. 



are evidence for template mismatch. Finally, we estimate the un- 
certainties from a set of 100 Monte Carlo simulations, where each 
spectrum has a an added perturbation consistent with random noise. 
Since we do not propagate the uncertainties of each spectral pixel in 
the data cube, the level to which a spectrum is perturbed is given by 
the robust standard deviation of the difference between the best fit 
and the original spectrum. During the Monte Carlo calculation the 
penalisation is turned off to produce realistic uncertainties. In the 
case of NGC524, typical derived errors are 10 km/s, 13 km/s, 0.03, 
0.04, 0.04 and 0.04 for V, a, h 3 , h 4 , h 5 and h 6 , respectively. Sim- 
ilarly, for NGC2549 typical errors are: 5 km/s, 7 km/s, 0.03, 0.03, 

0.04 and 0.04 for V, a, /13, hi, /is and hg, respecti vely. 

The SAURON kinematics was presented in lEmsellem et alj 

(2004). Here we use an extraction which includes up to he Gauss- 
Hermite moments and was do ne with Miles stellar templates 
JSanchez-Blazquez et alj2006t) in lCappellari et alj d2007t) . 

Figure [3] presents two-dimensional kinematics extracted from 
our NIFS observations. The nuclei of both galaxies show significant 
rotation, with both rotation velocity and velocity dispersion being 
higher in the more massive NGC524. In both cases, h$ and V are 
strongly anticorrelated, while the /14 maps are different: NGC524 
has a clear dip in the central 0'.'5, while the /14 map of NGC2549 is 
essentially constant and slightly positive over the whole field. 

Since the NIFS kinematics probe a completely different spec- 
tral range than the SAURON observations, it is necessary to ver- 
ify that these two data sets agree. The two sets of kinematics are 
widely separated in wavelength, and could therefore be expected 
to trace slightly different stellar populations with correspondingly 
different kinematics. Such population-dependent kinematics are not 
taken into account by our modelling process. Moreover, the abso- 
lute calibration of the velocity dispersion is sensitive to how well the 
intrinsic broadening of the template libraries (which for both data 
sets come from other instruments) mimics the instrumental disper- 
sion. To obtain a reliable model with a meaningful fit, the two data 
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Table 2. PSF of NIFS observations for NGC524 and NGC2549. 



galaxy 
image 



NGC524 
HST MGE 



NGC2549 
HST MGE 



fwhmjv 
fwhmg 
inten jy 



0.23 ±0.1 0.16 ±0.1 
1.25 ±0.1 1.36 ±0.14 
0.39 ±0.08 0.33 ±0.09 



0.17 ±0.03 0.20 ±0.02 
0.84 ±0.08 0.81 ±0.05 
0.53 ±0.08 0.47 ±0.05 



Notes - Table presents parameters of two-Gaussian fits to HST and 
deconvolved MGE model images. In the first row there are values for 
FWHM (in arcseconds) of the narrow Gaussian, in the second row there are 
values for FWHM (in arcseconds) of the broad Gaussian and in the third 
row is the intensity of the narrow Gaussian, where the sum of the intensities 
of the two Gaussians is one. 

sets must be mutually consistent within the assumed PSFs. In Fig.[4] 
we show a comparison between the velocity dispersions averaged 
on circular rings using bisquare weighting. The overlap is minimal; 
in radial direction NIFS data cover only about 2.5 SAURON pixels. 
Given the difference in the observational set ups and the uncertain- 
ties in the extraction of kinematics from such different spectral re- 
gions, it is remarkable to see such an agreement between the data 
sets. In the case of NGC524 we lowered the NIFS data for 4% to 
obtain the match with the SAURON data, while this is, in princi- 
ple, not necessary for NGC2549 NIFS data (a possible shift of 1% 
would still be consistent with the uncertainties). The agreement be- 
tween the velocity dispersion measurements is crucial for a robust 
determination of the mass of the black hole and, hence to constrain 
dynamical models of NGC524 we use the NIFS velocity disper- 
sions systematically lowered for 4% as shown in Fig. [4] 



4 DETERMINATION OF THE POINT SPREAD 
FUNCTION 

A robust estimate of the spatial resolution of our data is important 
for the construction of dynamical models, both in the sense of com- 
bining the data sets and in order to determine to which scales we can 
probe the internal dynamics of galaxies. Although the SAURON 
data were observed under natural seeing and the NIFS data were 
observed with LGS AO (see Section |2T2l without any PSF reference 
stars in the field, we use the same general method in both cases to 
determine the actual PSF. The FWHM of the PSFs of the SAURON 
observ ations for NGC524 an d NGC2549 are l'.'4 and l'.'7, respec- 
tively jEmsellem et all 2004), while NIFS PSFs derived below are 
listed in Table[2] 

Stars observed before and after the science frames show that 
there is minimal degradation of the AO correction during science 
observations (Section [2~2l which could arise from incorrect adapta- 
tion of the system to changes in focus of the LGS. This, however, 
does not guarantee that the AO correction of science frames (and the 
final resolution of the galaxy nuclei) will be of sufficient quality. 
Moreover, an accurate description of the PSF is required to com- 
pare the observations with the model predictions. Since there are 
no point sources in the field-of-view (stars or AGN sources such as 
nuclei or jet knots), the only way to estimate the achieved PSF is to 
compare the reconstructed NIFS images with images of higher res- 
olution: those obtained with the HST (for a list of various methods 
to estimate the PSF see |Paviesll2007t) . We first convolve the HST 
image with a concentric and circular double Gaussian description 
of the PSF, parameterised as the dispersions of the two components 
and their relative weight. The convolved image is then rebinned to 





1 1 1 1 


i i 1 

o NIFS _ 






HST 




i 

[ l 


MGE 




1 


HST convolved _ 




1 


MGE convolved^ 





















R [orcsec] 
NGC2549 



1 1 1 1 1 1 


1 1 1 

O NIFS 




— HST 




j — MGE 




— HST convolved 




.1 — MGE convolved . 

i 1 
t 1 

'A' 

I V 











R [orcsec] 

Figure 5. Comparison of light profiles along the major axis between NIFS 
and HST images for NGC524 (top) and NGC2549 (bottom). In both images 
NIFS profiles are shown with open symbols, HST profiles with black solid 
line, deconvolved MGE model profiles with red dashed line. HST and MGE 
profiles are convolved with appropriate PSFs from Table[2] Convolved HST 
and MGE profiles are presented by blue solid and orange dot-dashed lines. 



the same pixel size as the NIFS observation. This image is com- 
pared with the reconstructed NIFS image, and the parameters of the 
PSF are varied until the best matching double-Gaussian is found. To 
obtain the final values one has to quadratically add the PSF of the 
HST image to the measured PSF. The double-Gaussian parameter- 
ization of the PSF is to some extent degenerate: different combina- 
tion of the Gaussian parameters can reproduce the similar combined 
profiles. We estimate the uncertainties to the Gaussian parameters 
applying the method in different ways: we vary the initial condi- 
tions, we keep the centre of the images fix and let it vary freely, 
and we change the size of the NIFS maps limiting the effective area 
used in the comparison. Th e same method was used previously for 
SAURON and OASIS dat a jEmsellem et~aT]|2004l ; iMcDermid etal] 
l2006l : IShapiro et alj2006l) . 

Implicit in this method is the assumption that the reference im- 
age (here coming from HST) has significantly better spatial resolu- 
tion than the unknown image, such that the convolution kernel dom- 
inates the resulting image quality. The intrinsic spatial sampling of 
the HST/WFPC2 PC chip images is 0'.'0455 and the FWHM of the 
PSF of F702W and F814W images being used is 0'.'07 and 0"08, 
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Figure 6. Integrated flux from the determined PSF (using the HST images) 
for NGC524 (dashed line) and NGC2549 (solid line). The integrated flux is 
normalised such that the integral of the PSF for each galaxy is 1 . 



respe ctively, estimated using the Tiny Tim software dKrist &Hoold 
l200ll) . The final NIFS data-cubes have a pixel scale of 0'.'05 and in 
ideal conditions one could expect an AO-corrected PSF of ~ O'.'l 
FWHM as seen in Fig. Q] These values are very similar to the 
HST/WFPC2 characteristics, and the assumption that the PSF of 
the reference image is negligible may not strictly hold. We continue 
following two approaches: one uses the HST/WFPC2 images di- 
rectly, while the other uses a deconvolved model of the HST images 
as a PSF reference. In th e latter case, we use the Multi-Gaussian Ex- 
pansion (MGE) m ethod dMonnet et alJl992l : lEmsellem et al J 19941: 
ICappellarill20fB) to parameterise the HST images and deconvolve 
them analytically by their respective PSFs. The MGE method and 
the same deconvolved models are also used in the construction of 
the mass model use as input to the dynamical models of our galaxies 
(see Section lBTTl for more details). 

The individual observations of NGC524 and NGC2549 differ 
in the quality of the AO correction. If one overplots isophotes of 
the reconstructed images, one can see that some images are more 
peaked in the centre while others have distorted isophotes elongated 
in one direction. These differences arise from the changing LGS AO 
correction in response to the natural conditions. While one would 
ideally like to select only the best frames, there is a trade off be- 
tween the achieved S/N ratio and the size of the PSF. Inspecting 
individual data cubes by eye, we selected those which had regu- 
lar and concentric isophotes and merged them as described in Sec- 
tion |3.2| We verified that adding more lower quality frames worsens 
the measured PSF without significantly improving the S/N ratio. 
The final PSFs of the data used in this work are listed in Table|2] 

Figure [5] compares the PSF of our final NIFS data cubes for 
NGC524 and NGC2549 with the HST/WFPC2 and deconvolved 
MGE model images. They clearly show that our NIFS data are of 
lower spatial resolution than HST/WFPC2 images and, in this case, 
the MGE models are not really necessary for the PSF estimates, 
since the comparison with the deconvolved MGE model images 
yields very similar results to when using HST images directly. This 
is expected if the resolution of the data is sufficiently lower than the 
HST resolution in which case its PSF can be neglected. We there- 
fore assume the PSF derived from the HST images directly, rather 
than introduce the added step of using the MGE model, and refer to 
these values for the rest of the paper. 



For both galaxies, the narrow Gaussian is somewhat larger 
than the narrow component found with the tuning star (FWHM of 
0'.'23 and O'.'l 7 for NGC524 and NGC2549, respectively), while 
the poorer tip-tilt correction is most apparent in the larger width of 
the broad component (FWHM of 1'.'25 and 0'.'81 for NGC524 and 
NGC2549, respectively). NGC524 clearly has a poorer correction, 
but also uncertainties on the PSF components are larger, which is 
likely a consequen ce of core like light profile. As pointed out by 
iDavies et al.l d2008h . however, the Strehl ratio in LGS observations 
is not very dependent on the tip-tilt correction, since the flux within 
the core of the PSF does not change significantly if tip-tilt is present 
or not. We show the cumulative encircled energy of the PSF of our 
observations for NGC524 and NGC2549 in Fig. [6] In both cases 
there is approximately 40-50% of total flux within 0'.'2. The broad 
component of the PSF seems largely determined by the steepness of 
the galaxy's light profile, with the steeper profile of NGC2549 re- 
sulting in an overall narrower PSF and 90% of the total flux within 
0'.'6. For NGC524 90% of flux is achieved only within about 1". 



5 DYNAMICAL MODELS 

In this section we construct orbit superposition models defined by 
three integrals of motion and based on Schwarzschild's method 
dSchwarzscMloll979|;lRichstone & Tremainelll988l : lRix et aljl997l: 
van der Ma rel et al.l Il998l) . Our models are axisymmetric and 
contain further developments which include fits to the two- 
dimensional stellar kinematics and are adapted for more gen- 
eral surface-brightness distributions. Due to their generality, ax- 
isymmetric three-integral models are standard for M. determina- 
tions and ex ploration of the nuclear orbital structure in nearby 
galax i es (e.g.bebhardt et alj2003l;IValluri et alj2005l : lNowak et all 
l2007l . l2008tlGultekin et alj|2009l) . For more general triaxial mod- 
els see Ide Lorenzi et alj d2007h : Ivan den Bosch et al.l d2008l) . The 
implementation and adopted setup used here is described in 
Cappe ilari et al.l d2006l) . It consists of four steps: (i) the stellar light 
is parameterised, deprojected and converted to a stellar potential 
assuming axisymmetric shape and a stellar M/L; (ii) a represen- 
tative (dithered) orbit library evenly sampling across the observ- 
able space (the three integrals of motion and covering the luminous 
mass of the galaxy) is constructed; (iii) orbits are projected onto 
the observable space (sky positions and LOSVD parameters) tak- 
ing into account the PSF and apertures (Voronoi bins); (iv) weights 
to each orbit are determin ed using a non-negative least-squares fit 
dLawson & Hansonl 19741) . which when co-added reproduce the ob- 
served stellar density and kinematics in each bin. 



5.1 The mass model 

We used HST/WFPC2 F814W and F702W, and ground-based 
MDM / and R band images for NGC524 and NGC2549, re- 
spectively, to parameterise the s urface brightness dist r ibutions of 
our g alaxies with MGE models (Emsellem et al. 1994; Cappeilari 
120021) . Our MGE models we re corrected for galactic extinction fol- 
lowing Schlegel et alj dl998l) . as given by NASA/IPAC Extragalac- 
tic Data base (NED) and converted to a surface den sity in solar units 
using the WFPC2 calibration from Dolphin (2000), while assuming 
for / and R band absolute magnitude for the Sun of 4.08 and 4.42 
mag, respectively dBinnev & Merrifiel3ll998h . 

The MGE model for NG C524 was previously given in Table 
Bl of ICappellari et al. I d2006|) . We list the parameters of the MGE 
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Figure 7. HST WFPC2/PC images of NGC524 (left) and NGC2549 (right) 
with overplotted MGE models. Contours are spaced in half magnitude steps. 



Table 3. Summary of relevant properties for NGC524 and NGC2549 



Galaxy property 


NGC524 


NGC2549 


source. 


Morphological type 


50+ (rs) 


S0°(r)sp 


1 


Nuker 7 


0.03 


0.67 


2 


Nuclear type 


core 


cusp 


2 


Effective radius [arcsec] 


51 


20 


3 


cr e [km/s] 


235 


145 


3 


Distance [Mpc] 


23.3 


12.3 


4 


Inclination [degrees] 


20 


90 


5 



Notes - S ources (if two r eferences presen t, they refer t o galaxies respec - 
tive ly): 1 -limsellem et ai]l2004l); 2 -|Lauer et aI.ll2007l) ; lRest et alJfeOOll); 
3 - lEmsellem et alj ^200% ; 4 - iTonrv et alj feOOll) ; 2 - ICappellari et alj 
120061) . 



model for NGC2549 in Table|AT]of Appendix [A] and plot the mod- 
els over the HST images in Fig. [7] The minute differences between 
the MGE model and the image of NGC524 originate in the dusty 
central region of this galaxy. The galaxy isophotes are nearly round, 
and the MGE model was constructed to have axial ratios of all 
Gaussians equal to 0.95. NGC2549 is, however, an edge on galaxy 
and with a pecul iar isophotal shape. While i t is usually classified as 
a normal SO (e.g. lde Vaucouleurs et aljl99lh . this object has several 
characteristics which indicate a possible bar: its isophotes change 
from being st rongly boxy (Q',' 7 — 4"5) to strongly disky (4 — 20" 
and beyond) dRest et alj|200ll) , an 'X' shape distribution of light in 
the bulge (giving rise to strongly boxy isophotes) visible on HST 
images, and evidence for a ring-like structure at larger radii (vis- 
ible on ground based colour images). The variation from strongly 
boxy to disky isophotes is difficult to reproduce with only positive 
Gaussians, and hence the largest discrepancies (a few percent) in 
our MGE model arise in that region (4 — 5"). 

We assume that both galaxies can be well reproduced by an 
axisymmetric shape. These galaxies do not show isophotal twists 
and their kinematics have clear disk-like rotation dKrainovic et al.1 
2008). Since a proper dynamical treatment of bars is beyond the 
scope of this paper, we keep our axisymmetric assumption as the 
next closest representation of NGC2549. The bias introduced by 
modelling a likely barred galaxy using a model with a static, ax- 
isymmetric potential has so far not been well explored. What is 
clear, however, is that as we probe systems with lower black hole 
masses, bars become a common feature of the general dynamical 
make-up of the hosts. The degree to which the large-scale dynami- 
cal structure can affect the central dynamics is difficult to estimate, 
and as we shall discuss in Section [6] characterising the large-scale 
galaxy can help constrain M.. In the current absence of a better 
model, however, we will pr oceed as others have done using a st atic, 
axisymmetric model (e.g. Onken et al. 20071: Isiopis et alj2009h . 

Deprojection of the surface brightness of an elli psoidal body 
is formally n on-unique for all but the edge-on case jRvbicki|[l987l : 
iFranxl 1 19881) . In practice, determination of the viewing angles 
throug h stellar dynamical models is also found to be de gener- 



ate iKrainovic et alj2005l : lvan den Bosch & van de Venll2008l) . It is 
only possible to obtain an estimate of the inclination for such ob- 
jects by applying further assumptio ns, such as of the shape of the 
velocity ellipsoid (Cappellari 2008), u sing the ob servationally mo- 
tivated positive anisotropy ( Cappellari e t al. 2007). The accuracy of 
M. determinations from stellar dynamics are dependent on these 
effects, since one must choose an inclination to fix the potential. 
For the two galaxies considered here, there are reliable indications 
of their true inclinations. NGC524 has a large dust disk, which, if 



assumed to be intrinsic ally circular, gives an inclination of i ~ 20° 
(Cappellari et al. 2006). NGC2549 is close to edge on, which is vi s- 
ible both from photometry and kinematics dKrajnovic et al . 2008). 

While there are difference in stellar populations be- 
tween NGC524 and NGC2549, within each galaxy, the vari- 
ation line-strength indi ces are mild across the SAURON field 
(Kunts chner et al.ll2006l) . Kuntschner et al. (2009, in prep.) derive 
the age (single stellar population equivalent) within Re/8 apertures 
to be old (> 12 Gyr) and intermediate (5.5 Gyr) in NGC524 and 
NGC2549, respectively. Since M/L ratio is mostly a function of the 
Hp line-strength index and, hence, predominantly the age of stellar 
population (see Fig. 16 in ICappellari et al]|2006T) . which changes 
mildly across the fields, we assume in our dynamical models a con- 
stant M/L. We also assume that the dark m atter contributes only 



with a sma l l frac ti on in the centra l regions ( Gerhard et al. 2001 
Rusin et al.1 120031: ICappellari et al. 2006; Koop mans et alj 12006 



Thom as et al.ll2007l ; iBolton et alj|2008l) and do not explicitly add 



it to the models. 

Finally, in this paper we use distances from Tonrv et aljfeoOlh . 
adjusted for the Cepheid zero-point of iFreedman et al.1 (2001) 
putting NGC524 and NGC2549 at 23.3 and 12.3 Mpc, respectively. 



5.2 Stellar dynamical models 

Our Schwarzschild models were simultaneously fitted to the 
SAURON and NIFS kinematics and mass distributions parame- 
terised by the MGE models. Since the models are axisymmetric 
and by construction point-(anti)symmetric, we symmetrised the 
kinematic maps using kinematic position angles (PAfci„ = 1° 
and 40° for NGC524 and NGC2549, respectively) determined in 
Capp ellari et ai]( l2007h . The symmetrisation uses the mirror-(anti)- 
symmetry of the kinematic maps, such that kinematic values from 
for positions ((x,y,), (x,-y), (-x,y), (-x,-y)) were averaged. Since the 
Voronoi bins have irregular shapes and they are not equally dis- 
tributed with respect to the symmetry axes, we average the four 
symmetric points, and, if for a given bin there are no bins on the 
symmetric positions, we interpolate the values on those positions 
and average them. We keep the original errors of these bins in or- 
der not to underestimate the parameter uncertainties. The assump- 
tion here is that the uncertainties of the kinematic parameters in 
bins at positions (x,y,), (x,-y), (-x,y), (-x,-y) are similar, which is 
reasonable given that the light falls off symmetrically in all four 
quadrants. We exclude SAURON bins in the central l'.'3 (9 bins), 
since these overlap with the high resolution NIFS observations. We 
constrain the dynamical models fitting the kinematics by up to 
Gauss-Hermite moments for both data sets. 
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Figure 8. Schwarzschild dynamical models and the determination of the best fitting parameters for NGC524 (left) and NGC2549 (right). Each symbol is a 
dynamical model specified by given M/L and M. . The agreement between the data and the models are described by overploted contours A\ 2 = X 2 ~ Xmi„ 
contours showing 1 , 2 and 3<r levels for two parameters. Further contours are spaced by a factor of 2. Open (red) circles mark the best fitting models. Top panels 
show grids of models constrained by both SAURON and NIFS kinematics, middle panels show grids of models constrained by SAURON only kinematics and 
bottom panels show models constrained only using NIFS data. 



Each dynamical model is defined by two free parameters: M/L 
and M. . The orbit library is constructed for a given M. at a given 
(expected) M/L and consists of 444528 orbits, which are bundled 
in groups of 6 3 = 216 before the linear orbital superposition. It is, 
however, not necessary to compute the orbit library for each pair 
of (M/L, M.). The predictions of the orbit library can be scaled to 
construct models for different M/L ratio s. We use a modest amoun t 
of regularisation A = 10 (as defined in lvan der Marel et al.lll998h . 
Figure [8] shows the main results of this study: for a set of (M/L, 
M.) pairs defining individual models overploted are contours of x 2 , 
showing the agreement between the kinematic data and the model 
predictions. 

In the case of NGC524, the best fitting model has 
M.=(8.3l f 7 3 ) x 10 8 M and M/L = 5.8 ± 0.4 (/ band). Note that 
the M/L estimate for the almost face-on NGC524 i s strongly depen- 
dent o n the assumed inclination (see Figure Al in ICappellari et alj 
2006). For NGC2549, the best fitting model is the one with 
M.=(1.4±?;§) x 10 7 M and M/L = 4.7 ± 0.2 (R band). The 
quoted uncertainties for M. are 3a values marginalised over the 
M/L and correspond to one degree-of-freedom (Ax 2 = 9). Due to 
numerical uncertainties inherent to Schwarzschild' s method we ad- 
vise using 3a errors only and we do not measure the uncertainties 
at 1 and 2 a level s (for a discussion on top ology of x 2 contours 
see Appendix A of Ivan der Marel et apfT998T) . Since errors are ex- 
pected to scale as yj Ax 2 JPress et al.l 1992h one can estimate 1 and 



2 a confidence levels as approximately 1/3 a nd 2/3 of the 3a values , 
respectively. Similar approach was used bv lGebhardtetal.U2003l) 
to estimate 1 a confidence levels of studies that only provide 3a. 
In addition to models fitted to both the SAURON and NIFS data, 
we also constructed models constrained with only NIFS and only 
SAURON data, in order to asses the relative importance of these 
data sets. This exercise demonstrates that the best-fitting parame- 
ters are consistent within 3a level between models constrained fit- 
ting different data sets. 

More importantly, this demonstrates that the SAURON data, 
whilst being relatively insensitive to the black hole mass due to its 
lower spatial resolution, provides a critical constraint on the M/L 
by virtue of its large-scale spatial coverage. The M/L and M. are 
not independent, and the drastic reduction in uncertainty of the M/L 
provides a significant reduction in the uncertainty of M.. Without 
the SAURON data, the uncertainty for M. doubles, while errors 
on M/L increase between 5 to 10 times. The NIFS data are how- 
ever crucial to determine M., especially for low mass M. when 
Rsph becomes much smaller than a SAURON resolution element; 
'SAURON only' model for NGC2549 can only give an upper limit 
to M. , and NIFS data are necessary to determine it robustly (see 
also lShapiro et alfcOQrj) . 

In Figs.l9land[T0l we show the residuals between the data and 
the best fit models for NGC524 and NGC2549, respectively. The 
residuals were calculated by subtracting the predictions of the best 
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Figure 9. Residuals between the data and the best fit Schwarzschild model for NGC524 presented by NIFS (top) and SAURON (bottom) data. The residuals 
were obtained subtracting the model from the data and dividing by the errors. From left to right the columns show: mean velocity V, velocity dispersion cr, and 
Gauss-Hermite moments /13 — /ig . Overplotted contours are surface brightness is ophotes. Circles i n the c enter of the SAURON maps show bins excluded from 
the fit. North is up and East to the left. SAURON kinematic maps can be found in lEmsellem et al.1 )*2 004l) . 
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Figure 10. Residuals between the data and the best fit Schwarzschild model for NGC2549 presented by NIFS (top) and SAURON (bottom) data. The residuals 
were obtained subtracting the model from the data and dividing by the errors. From left to right the columns show: mean velocity V, velocity dispersion cr, and 
Gauss-Hermite moments — h^. Overplotted contours are surface brightness isoph otes. Circles in the center of the SAURON maps show bins excluded from 
the fit. North is up and East to the left. SAURON kinematic maps can be found in Emsel lem et al 



fit model from all kinematic constraints (V, a, hi - ha) and dividing 
by the associated errors. The stretch in colour is chosen to show 
the goodness of fit relative to the kinematic errors: the extremes 
bound 3 times the uncertainty. Most of bins, for both NGC524 and 
NGC2549, have green, yellow or light blue colours which mean 
that the residuals are within the kinematic errors. Somewhat higher 
residuals (> 2a level) are systematically visible only for higher 
Gauss-Hermite moments (h$, he)- 

The general nature of the three integral models allows us to 
investigate the orbital anisotropy of our best fit models. We do this 
by plotting different coordinates of the velocity dispersion tensor. 
We first define the ratio of radial to tangential velocity dispersions, 
oy/<7 f in Fig. QT] where of = (ag + a^,)/2, in spherical coor- 



dinates (r, 8 and (f>). Since includes only random motion, an 
isotropic system will have a r /a t = 1. This measure of anisotropy 
has been used in some previous studies of n uclear orbital structure 
dGebhardt et al l 120031 ; Ishapiro et alj [20061 : iHoughton et"ai] 120061 : 
Giilte kin et alj 120090 . but we note that it is mostly suited for de- 
scription of spherically symmetric objects. NGC2549 is an edge on 
galaxy and clearly not spherical, while the low inclination and the 
projected flattening of MGE Gaussians determin e the intrinsic flat - 
tening of NGC524 to be q, ~ 0.41 (eq. (9) of ICappellarill2002l) . 
and again not a spherical galax}0. For these reasons we also make 



5 The regions under the direct influence of SMBHs tend to be spherically 
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use of the components of the velocity dispersion tensor in cylindri- 
cal coordinates (R, cj> and z): or, 0$ and o z . We define anisotropy 
parameters /3 X — 1 — (o z /or) 2 and 7 = 1 — (o^/or) 2 . The 
first parameter, /3 Z , describes the shape of the velocity ellipsoid in 
(vr, v z ) plane, a plane including the symmetry axis plane of an 
axisymmetric galaxy, where j3 z = if the cross-section of the ve- 
locity ellipsoid is a circle. The second parameter, 7, describes the 
shape of the velocity ellipsoid in (vr, v$), a plane orthogonal to 
v z , which can be identified with the equatorial plane in an axisym- 
metric galaxy. Again, if 7 = 0, the shape of the velocity ellipsoid 
in this plane is a circle. Departures from the circular shape of the 
velocity ellipsoid quantify the bias towards radial ( f3 z > 0, 7 > 0) 
or tangentially ( f3 z < 0, 7 < 0) biased orbital distributions. 

In Fig. QT] we plot these three anisotropy parameters for 
NGC524 and NGC2549 as as function of radius and position angle 
within each galaxy. The galaxies show similar orbital structures at 
large radii, but markedly different orbital structure close to the black 
holes. In the spherical coordinates, orbits in NGC524 are clearly 
tangentially anisotropic. The significant drop in ay/o> occurs at 
the radius similar to the R ap h, which also corresponds well to the 
radius within which there is a significant drop in measured values 
of the ti4 Gauss-Hermite moments. To see this effect in cylindrical 
coordinates one has to keep in mind the definition of the coordinate 
systems, where the r and R coordinates have the same orientation 
on the major axis, but a perpendicular on the minor axis, hence the 
difference between the minor and major axis in {3 Z parameters. The 
equatorial plane view of 7 anisotropy values are also consistent whit 
this picture. The orbital bias is tangential, even without the stream- 
ing motion present in anisotropy parameters. Note that for the minor 
axis 7=0 by symmetry for an axisymmetric galaxy. For NGC524 we 
can conclude that both spherical and cylindrical coordinate systems 
give a clear picture of orbits being mostly tangentially biased in the 
central C/.'5. At larger radii, where the flatness of the galaxy is more 
dominant, the anisotropy parameters defined in the cylindrical co- 
ordinates shows that o z < or, but or ~ 0$, which means that the 
velocity ellipsoid in N GC524 is very close to oblate (as shown in 
ICappellari etai]|2007l) . 

NGC2549 shows a somewhat different picture in the nucleus. 
The spread of anisotropy in the cr r /crt is inconclusive regardless 
the anisotropy (or possibly weakly tangentially anisotropic). On the 
other hand, f3 x clearly shows that anisotropy is very similar along 
both major and minor axis, and in generally o z < or, without a 
notable change towards the central black hole, as it is visible for 
NGC524. A possibly important difference is that for NGC524 the 
resolution of our observations is several times less than the K sp h, 
while in NGC2549 they are only comparable, and a possible tan- 
gential anisotropy structure around the SMBH is not visible. In the 
equatorial plane, 7 values along most position angles are basically 
consistent with zero. The values of /3 Z and 7 imply that NGC2549 
has oblate velocity ellipsoid. The anisotropy of NGC2549 over the 
SAURON field of view is in an excellent agreement with the re- 
sults from Jeans m odelling with con stant anisotropy (/3 Z = 0.17 
and 7 = 0, Fig. 5 in Ca ppeUarill2008h . 



symmetric allowing a representation of the moments of the velocity ellipsoid 
in spherical coordinates. 
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Figure 11. Radial profiles of anistropy of Schwarzschild dynamical models 
for NGC524 (left) an NGC2549 (right). Top: anisotropy parameterised in 
spherical coordinate system by ay /<rt . Middle: anisotropy parameterised in 
cylindrical coordinate system by f3 z . Bottom: anisotropy parameterised in 
cylindrical coordinate system by 7 (see text for definitions). Anisotropy was 
measured at different polar angels defined in the meridional plane, ranging 
from 15° (close to the equatorial plane, red line) to 75° (close to the sym- 
metry axis, blue line). The change in colours shows the increase in the angle 
from the major to minor axis (red to blue colours). Dashed vertical lines 
show cr of the narrow component of the LGS PSF and dot-dashed vertical 
lines show R sp h using the derived M. . 



6 DISCUSSION 

6.1 Relative importance of NIFS and SAURON data 

The comparison between the grids of models constrained by NIFS 
and SAURON, only NIFS and only SAURON data (Fig.[§J show the 
importance of combined large scale and high resolutions data sets. 
Stellar dynamical models constrained by high-resolution data with 
small field-of-view are not sufficient to robustly determine M.. The 
main reason is that they are not able to constrain the orbital distri- 
bution and hence the total M/L of the galaxy. This reflects the need 
to establish the importance of stars on radial orbits, which may pass 
close to the SMBH, whilst spending most of their time at much 
larger radii. In the lower panels of Fig.[8]the models constrained by 
NIFS data alone show a clear degeneracy between the M/L and M. 
where low M/L is increasing M. and vice versa. While the formal 
best fitting M. are not very different from the one determined using 
both NIFS and SAURON data, the uncertainties are much larger: 
a factor of a few in M. to about an order of magnitude increase 
in M/L uncertainties. The SAURON data (middle panel on Fig. [8} 
constrain the M/L of our galaxies, becaus e they cover mo r e than 
1 R e , which is consistent with finding of iKrainovic et al. I d2005h 
that large scale data carrying information on the orbital structure 
over approximately l_R e are necessary to constrain the distribution 
function. The implications are clear: to accurately measure the cen- 
tral orbital structure and gravitational potential (i.e. including M.) 
in galaxy nuclei requires IFU kinematics on both small and large 
scales. Similar results, using high resolution and large scale IFU 
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data, were reported in ICopin et alj J2004I) and lShapiro etafl J20061) . 
except that in the case of NGC2549, a smaller galaxy compared to 
the galaxies in these studies, the high resolution data are actually 
crucial to determine the lower limit to M.. 

The sphere of influence for NGC524 is R sp h ~ 0'.'6 and for 
NGC2549 R sp h ~ 0'.'05 using the masses derived in this study. In 
the case of NGC524 the R sp h is well resolved by the LGS observa- 
tions, while the SAURON data have the sigma of the Gaussian PSF 
quite similar (0'.'6). In the case of NGC2549, the sigma of the nar- 
row Gaussian of the LGS observations (0'.'07) is similarly just larger 
than the R sph . While neither SAURON data for NGC524 nor NIFS 
data for NGC2549 formally resolve respective R sp h, they do con- 
strain the M. . The reason likely lies in the combination of high S/N 
of data, IFU coverage of the nuclear region and a significant part 
of the galaxy and that sigmas of the PSFs are similar to the R sp h, 
partially resolving the spheres of influence. This finding is impor- 
tant for future dynamical studies of SMBHs as well as for studies 
of M. — <r e relation. A discussion of the biasses introduced when 
M. measurements based on data which do not fully resolve R 3p h 
are exclud ed from determ i ning t he properties of M. — <r e relation, 
is given in iGultekin et alj |2009). We, however, note that likely the 
crucial contribution for determination of M. in marginally resolved 
cases comes from two-dimensional coverage of the kinematic con- 
straints. 

In order to illustrate that our models can robustly measure 
M. in these galaxies, in spite of marginally resolving R sp h, in 
Fig. [12] we show clearly visible difference between the data and 
various models. Here we concentrate on the models constrained by 
SAURON and NIFS data for NGC2549 and only SAURON data for 
NGC524. In the latter case we want to stress that even SAURON 
data, without the help of high spatial resolution data are able to 
distinguish between models with too large and too small M. . In 
the upper panels of the Fig. [T2] we show the difference between 
velocity dispersion maps for three models (SAURON + NIFS con- 
straints) at the best fitting M/L for different M. in NGC2549: the 
best fitting model, a model with too small and a model with too 
large M.. The difference between models is mostly evident in the 
central a peak: a too low M. underpredicts, while a too high M. 
overpredicts it. The contribution of the NIFS data in constraining 
the lower limit on the M. in NGC2549 is not only statistical and 
visible in \ 2 values, but can also be judged by eye. In the lower 
panels of Fig.[l2]we plot a comparison between the velocity disper- 
sion maps for data and three models (SAURON only constraints): 
the best fitting model, a model with too small and a model with 
too large M. . Again, the only relevant change between the model 
maps occurs in the central few pixels, with the largest change in the 
central pixel, confirming that even using only SAURON data can 
constrain the M. in NGC524, although the uncertainties are larger 
than when NIFS data are added. 

6.2 Accuracy of M. determination on the PSF 

Observing nearby galaxies with AO systems is not practical unless 
it is possible to use the LGS without relaying on nearby tip-tilt stars, 
such as in this work. The correction of the PSF, however, in simi- 
lar cases will depend on the properties of the central light profiles. 
From that point of view we want to establish to what accuracy it is 
possible to measure M. covering a range of corrections achievable 
with current LGS systems. 

We tested this by constructing two simple model galaxies de- 
scribed by a de Vaucouleurs R 1 ' 4 law, of total mass 5 x 10 10 M Q 
and size R e = 27", also adding an SMBH of 2% and 1% of the 



galaxy mass in each of them. When put at a distance of Virgo (16.5 
Mpc), the SMBH in these model galaxies had estimated R sp h of 
0'.'27 and 0'.'13, respectively. We assumed a sph erical symmetry a nd 
used Jeans models (the JAM implementation of lCappel lari 2008) to 
predict the second kinematic moments as a function of radius, given 
the light of our model galaxies parametrized with the MGE method 
and a realistic PSF. We made a 100 different model predictions of 
second moments (Vrms), each with a different PSF. The PSF was 
parameterised as a double Gaussian where the FWHM of the broad 
component is kept fixed to a (natural) seeing of 0'.'8, while the nar- 
row component model was varied from 0'.'04 (approximating the 
diffraction limit at 8-10 meter telescopes in the near-infrared) to 
0'.'5. The relative flux ratio of the components was such that their 
sum is equal to unity and we varied the flux of the narrow compo- 
nent from 0. 1 to 0.5. 

We then added a constant error of 10 km/s to the derived V rmis 
and fitted them with the spherical Jeans models, but now varying the 
M. (100 different values centred on the true value) and account- 
ing for M/L by scaling the predicted Vrms to the model. After the 
best fit was found we determined the range of M. within 3<j level 
(Ax 2 = 9). Figure [T3l illustrates this process. The symbols repre- 
sent the Vrms of the galaxy model which are fit with Jeans models. 
The model with the best fit M. is shown with a dashed line, while 
the two lines above and below the best fit model are the models that 
are at A\ 2 = 9 level away from the best fit M.. For each PSF 
we calculate the difference in values of M, for these two models. 
This number, expressed at a fraction of the input M., is an estimate 
of the accuracy to which the M. can be determined given the PSF 
properties. 

The results for all PSFs are shown in Fig. [14] As expected, 
regardless of the size of the SMBH, the best accuracy is achieved 
for small FWHM and high flux of the narrow-Gaussian component. 
The size of the SMBH, however, is important for the absolute ac- 
curacy: bigger SMBHs will be determined with higher accuracy. 
A typical PSF achievable with the LGS systems, such as in this 
study, of 0'.'2 FWHM and 0.4 intensity of the narrow component of 
the PSF, yields an uncertainty of about 40% for R sp h=0'.'27 and an 
uncertainty of about 70% for R sp f 1 =0'.'l3. The uncertainties of our 
determinations for NGC524 and NGC2549 are approximately 34% 
and 54%, respectively, which are in a good agreement with the pre- 
dicted uncertainties of 40% and 55% from these simple simulations 
(diamonds in Fig. 1141 . 

Real galaxies are neither pure R 1/4 laws nor spherically sym- 
metric objects; this exercise is a simplistic in many ways, however, 
it offers an insight to what improvement one can expect from the 
current LGS AO observations for SMBH studies. In a very good 
case of a large SMBH, when the PSF correction is approaching 
diffraction limit and the highest possible Strehl, it can be expected 
that the overall uncertainty approaches 25-30%. On the other hand, 
it is highly encouraging that even masses of small SMBHs can ro- 
bustly determined for quite low PSF corrections of e.g. 0'.'4 and 
Strehl of ~ 0.2. In general, it can be expected that for a decent AO 
correction, a typical error in M. in a nearby galaxy will be of the 
order of a factor of 2. 

This simulations suggest that there is an 'instrumental' limit 
on individual measurements of 30% which will translate into the 
intrinsic scatter of M. — <r e relation. Note that due to simple as- 
sumptions involved, the simulations can only be used as a guideline 
for the expected relative improvement in M. accuracy when using 
AO corrections, and that the percentages given in these simulations 
are dependent on the assumed kinematic errors. It, however, does 
suggest that high signal-to-noise integral field data obtained with 
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Figure 12. Comparison between data and model velocity dispersion maps for (top) NGC2549 and NGC524 (bottom). From left to right: cr (symmetrised) 
and model prediction maps for different M. (at the best fitting M/L). The best fitting model is shown in the map adjacent to the observed a map, followed by 
a model with a too small M. and a model with too large M. , respectively. The models shown here for NGC2549 were constrained with both S AURON and 
NIFS data, while the models shown here for NGC524 were constrained using only SAURON data. Note the colour change in the central few pixies which are 
a consequence of change in M. . North is up and east to the left. 



LGS AO are capable of robustly measuring M. in nearby galaxies 
and, hence, determining the extent (low masses), the shape, as well 
as increasing limit on the accuracy of the intrinsic scatter of M. — a 
relation. 



6.3 Nuclear orbital structure and formation scenarios 

Whil e both NGC524 ad NGC2549 are classified as fast rota- 
tors l lEmsellem et ai] |2007h , and seem to have a similar shape but 
viewed at different inclinations, there are some significant differ- 
ence between them. NGC524 is more massive and has luminosity 
weighted stellar populations older than NGC2549. The light profile 
of NGC524 shows a shallow core at radii smaller than (/.'3, and of 
NGC2549 a cusp at the HST resolution. Masses of the black holes 
are different, and our models suggest the orbital structures around 
them are also different. These characteristics point out to a different 
evolution of these objects. 

It has been suggested that the existence of cores in light pro- 
files of massive galaxies is connected to the formation process 
of nuclei , where high nucl e ar densities are dest royed in gas poor 
mergers dFaber et al]ll997l ; iHopkins et alj|2009l) . The main point 
here is that during collisions of galaxies without significant amount 
of gas, the process of black hole merging removes a significant 
number of stars with orbits intersecting the newly formed black 
hole binary. The binary shrinks ejecting stars mostly on radial or- 
bits that approach the black holes and only stars on tangential or- 
bits remain in a depleted nucleus dMilosavlievic & MerrittllToOll ; 
IMilosavlievic et al.ll2002h This process goes in several steps where 
both classical and relativistic effects have to be taken into account 



(e.g. IMilosavlievic & Merrittll2003l ; iMerritt et~al]|2007l) A different 
scenario suggests that mergers of galaxies with significant amount 
of gas lead to a creation of cuspy nuclear ligh t profiles, which are 
created as a consequence o f a nuclear starburst dMihos & Hernquistl 
1 1994 IHopkins et 311 2009 ) . To some extent this latter scenario in- 
cludes accretion of cold intergalactic gas. These processes are sum- 
marised and discussed in more details in Ko rmendv et al. (2008) 
who put forward a hypothesis that "the last major merger that made 
core ellipticals was gas poor, whereas the last major merger that 
made coreless ellipticals was gas rich and included a substantial 
starburst". 

In all these respects the growth of black holes is related to the 
formation of the nuclei in early-type galaxies. It is somewhat un- 
fortunate that both growth through binary mergers and accretion of 
gas leave similar imprints in the orbital structure: a tangentially bi- 
ased anisotropy of the velocity ellipsoid. The difference between 
the scenarios, however, could be found in the specific angular mo- 
mentum content of nuclei, if gas poor merges are also responsible 
for produci ng slow rotators, whi le gas rich mergers create fast ro- 
tators (e.g. lEmsellem et alj |2007). If this is true, than slow rotators 
with cores should have tangential anisotropy in their nuclei. In that 
case, products of gas poor mergers are slow rotating galaxies, with 
light deficits and tangential anisotropy in the nuclear orbital distri- 
butions (within the core). Given that NGC524 is a fast rotator with 
these characteristics, it is clear that there is no a simple formation 
scenario. Finding, however, a slow rotator without a core, or tangen- 
tially biased orbits, will argue against gas poor merger scenario as 
outlined above. Given that kinematically decoupled cores with rel- 
atively higher specific angular momenta are often present in slow 
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Figure 13. Two examples of the simulations used to test the influence of the 
PSF on the accuracy of the M. estimates. Top: Predictions for model with 
a poor PSF given by FWHM and the flux of the narrow component of C'45 
and 0.17, respectively. Botom: Predictions for model with an excellent PSF 
given by FWHM and the flux of the narrow component of C/.'14 and 0.63, 
respectively. In both panels, the symbols represent the V rms of the galaxy 
model which are fit with Jeans models. The best fit model is given by the 
dashed line, and the two lines above and below the best fit model are models 
that are at A\ 2 = 9 level away from the best fit M, 



rotators ( Kraj novic etafl l2008), it would be interesting to see the 
orbital structure in the vicinity of their central black holes. 

Within present scenarios, the fact that NGC524 has a core and 
very old stellar population, but also a prominent dusty disk, and 
has evidence of being a spatially constrained post-starburst galaxy 
(Kuntschner et al. in prep, Shapiro et al. in prep), suggest that in 
this galaxy first a core was created (i.e through a binary black hole 
merger), but soon after that accretion of gas created the disk and 
specified the angular momentum of the galaxy, while also possi- 
bly contributing to the final growth of the black hole, keeping or 
enhancing the nuclear tangential anisotropy. 



7 CONCLUSIONS 

In this work we presented observations of two early-type galax- 
ies (NGC524 and NGC2549) with laser guide star adaptive optics 
obtained at GEMINI North telescope using the NIFS integral field 
spectrograph. The purpose of these observations was to measure the 
kinematics within the nuclei of these galaxies and, in combination 
with previously obtained wide field observations with SAURON 
at William Herschel Telescope, to determine the masses of central 



Figure 14. Uncertainty of M. determination relative to the flux and size 
(FWHM) of the narrow component of the PSF in the case of a model galaxy 
with SMBH whose R sph ~ 0'.'13 (top), and a model galaxy with SMBH 
whose R sp h ~ 0''27 (bottom). Solid points show the parameters of the 
narrow-Gaussian component of the PSF used in construction of Jeans mod- 
els. Over-plotted lines are contours of the equal uncertainty to which M, 
can be determined using the given PSF. The given numbers are the percent- 
age of the true value of the M. spaced in 4.5%. Filled diamonds mark the 
PSF values for our observations of NGC2549 (top panel) and NGC524 (bot- 
tom panel) for comparison. 



black holes by means of stellar dynamical modeling. A speciality of 
these observations was that there were no natural guide stars for the 
selected objects and they were performed using an 'open loop' fo- 
cus model. In this mode, the nucleus of the galaxy is used as a sub- 
stitute for the natural guide star, but since it is not bright or distinct 
enough, the variations of the laser guide star focus are not tracked in 
real time, but varied following a geometric model. We selected two 
galaxies with different light profiles in order to test the dependence 
of the AO corrections on the cusp/core nature of the nucleus. 

We estimated the seeing correction of the the LGS system by 
finding the PSF that, when used to convolve an HST/WFPC2 im- 
age or deconvolved MGE model, best matched the reconstructed 
NIFS images. The deconvolved MGE models were useful to verify 
that the derived PSF estimates were reliable, since the LGS AO ob- 
servations approach the spatial resolution of HST, making it poten- 
tially difficult to neglect the intrinsic HST image quality. We param- 
eterise the PSF with a double-Gaussian, consisting of a narrow and 
broad component. The final full-width-half-maximum PSFs of our 
observations are Of! 23 and l'.'25 for NGC524 and 0"17 and 0'.'84 
for NGC2549, with the narrow component flux contribution of 0.39 
and 0.53 for NGC524 and NGC2549, respectively. These values 
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demonstrate that the laser correction works well, giving a signif- 
icant narrow component close to diffraction limit, to some extent 
regardless of the central light profile. The broad component of the 
PSF seems to be also dependent on the steepness of the light profile. 
This can be understood in terms of how much low-order (tip-tilt) 
correction can be made using the nucleus, which directly depends 
on the total brightness and the steepness of the nuclear light profile. 
In addition, the uncertainties in the estimate of the PSF parame- 
ters depend on the steepness of the profile. Our observations con- 
firm that the 'open loop' focus model provides a robust and stable 
method by which a range of objects with faint, extended light pro- 
files can be successfully used as tip-tilt references for LGS AO. The 
steeper light profiles, however, will yield a larger contrast within the 
field-of-view of the wave front sensor, and result in the overall bet- 
ter PSF, making the 'open loop' method more suitable for smaller 
early-type galaxies. 

Our observations were taken in K band, and we use the CO 
band head to derive the stellar kinematics. We compare these kine- 
matics with the optical and large scale SAURON kinematics and 
find that they match well. Both data sets are used to constrain three- 
integral Schwarzschild dynamical models. The models are axisym- 
metric, numerical and based on an orbit superposition method. For 
NGC524 they yield M.=(8.3±?;£) x 10 8 M and for NGC2549 
M.=(l-4i?;a) x 10 7 M . We showed that large-field IFU data 
are crucial for constraining the global M/L and contribute to the 
determination of the M. in cases when black hole's sphere of influ- 
ence is of similar scales as the obtained spatial resolution. If only 
the high-resolution data of small field-of-view are used to constrain 
the models, the best fitting values are similar, but the uncertainties 
are considerably larger: up to 2 times for M. and 5 times for M/L. 
The combination of both large scale (SAURON) and high resolu- 
tion (NIFS) data is therefore essential to provide useful constraints 
on M.. 

The nuclear orbital distribution of the best fitting models is dif- 
ferent for these two galaxies. NGC524 has the clear signature of a 
tangentially anisotropic velocity ellipsoid in the central f/.'5, which 
corresponds to the SMBH sphere of influence, the change in the 
light profile towards a central core and change in the shape of the 
LOSVD. In the case of NGC2549 more significant is that o z < or, 
and or ~ over the whole field-of-view covered with the kine- 
matic data. The difference in the nuclear dynamics between these 
two galaxies suggest that they had different formation processes, 
likely linked with the growth of the central black holes. The lack of 
tangential anisotropy next to the SMBH in NGC2549 could, how- 
ever, be a consequence of barely resolving its sphere of influence. 

We constructed two spherical Jeans models of galaxies with 
SMBHs of different masses and different spheres of influence and 
explored a relative improvement one should expect in determina- 
tions of M. in these galaxies given a range in seeing parameters 
typical for present LGS AO system. We find that for PSF approach- 
ing diffraction limit and high Strehl ratios the uncertainty in M. 
approaches 25-30%, but for a more typical PSF it is about 50%. 
Our results suggest that systematic studies of SMBHs in the nearby 
galaxies with LGS AO are able to decrease the intrinsic scatter and 
determine the extent and shape of M. — <r e relation. 
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Table Al. Parameters of MGE models for NGC2549 
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